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Abstract 

Recently, an analytic method was developed to study in the large limit non-hermitean 
random matrices that are drawn from a large class of circularly symmetric non-Gaussian 
probability distributions, thus extending the existing Gaussian non-hermitean literature. One 
obtains an explicit algebraic equation for the integrated density of eigenvalues from which 
the Green's function and averaged density of eigenvalues could be calculated in a simple 
manner. Thus, that formalism may be thought of as the non-hermitean analog of the method 
due to Brezin, Itzykson, Parisi and Zuber for analyzing hermitean non-Gaussian random 
matrices. A somewhat surprising result is the so called "Single Ring" theorem, namely, that 
the domain of the eigenvalue distribution in the complex plane is either a disk or an annulus. 
In this paper we extend previous results and provide simple new explicit expressions for 
the radii of the eigenvalue distiobution and for the value of the eigenvalue density at the 
edges of the eigenvalue distribution of the non-hermitean matrix in terms of moments of the 
eigenvalue distribution of the associated hermitean matrix. We then present several numerical 
verifications of the previously obtained analytic results for the quartic ensemble and its phase 
transition from a disk shaped eigenvalue distribution to an annular distribution. Finally, we 
demonstrate numerically the "Single Ring" theorem for the sextic potential, namely, the 
potential of lowest degree for which the "Single Ring" theorem has non-trivial consequences. 
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1 Introduction 



There has been considerable interest in random non-hermitean matrices in recent 
years. Possible applications range over several areas of physics[0, ^ ^. For some 
recent reviews see . One difficulty is that the eigenvalues of non-hermitean matrices 
invade the complex plane, and consequently, various methods developed over the years 
to deal with random hermitean matrices are no longer applicable, as these methods 
typically all involve exploiting the powerful constraints of analytic function theory. 
(See in particular the paper by Brezin, Itzykson, Parisi, and Zuber [^.) In two 
of us proposed a "method of hermitization" , whereby a problem involving random 
non-hermitean matrices can be reduced to a problem involving random hermitean 
matrices, to which various standard methods (such as the diagrammatic method]^, 
or the "renormalization group" method [H, ^ |lll) can be applied. An idea similar 
to the "method of hermitization" was expressed independently in 0. 

To our knowledge, the literature on random non-hermitean matrices [|], ^ has 
focussed exclusively on Gaussian randomness. For instance, it has been known for 
over thirty years, from the work of Ginibre[0, that for the Gaussian probability 
distribution P{4>) = {1/Z)exp (— A^tr0"'^0) (here, as in the rest of this paper, denotes 
an X complex random matrix with the limit N —>■ oo understood), the density of 
eigenvalues of is uniformly distributed over a disk of radius 1 in the complex plane. 

Analytic determination of the density of eigenvalues of a non-Gaussian probability 
distribution of the form 

p(0) = 1 -A^try(0t<^) ^ 
Zj 

where V is an arbitrary polynomial of its argument, was given for the first time in 
IQ. Based on the method of hermitization, it was shown in that by a simple trick, 
the desired density of eigenvalues could be obtained with a minimal amount of work, 
by judiciously exploiting the existing literature on random hermitean matrices. 

Due to the symmetry of -P(0) under the transformation — e*"0, the density 
of eigenvalues is obviously rotational invariant. It was shown in |^ that the class of 
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probability distributions of the form ( |1 . 1| ) exhibits a universal behavior in the sense 
that whatever the polynomial V was, the shape of the eigenvalue distribution in the 
complex plane was always either a disk or an annulus. This result was referred to in 
m as the "Single Ring Theorem" . 

In a certain sense, the formalism developed in ^ may be thought of as the analog 
of the work of Brezin et al. for random hermitean matrices ; they showed how the 
density of eigenvalues of hermitean matrices (p taken from the probability distribution 
P((/?) = {l/Z)exp[~NtTV{ip)] with V an arbitrary polynomial can be determined, 
and not just for the Gaussian case studied by Wigner and others [|l3l, in which V = 



(l/2)tr(y9^. An important simplifying feature of the analysis in is that P{ip) depends 
only on the eigenvalues of ip, and not on the unitary matrix that diagonalizes it. 
In contrast, the probability distribution (|l.lj ) for non-hermitean matrices depends 
explicitly on the GL{N) matrix S used to diagonalize = S^^AS, and S does 
not decouple. Remarkably however, for the Gaussian P{(f)), Ginibre [jl2| managed 



to integrate over 5* explicitly and derived an explicit expression for the probability 
distribution of the eigenvalues of 0. Unfortunately, it is not clear how to integrate 
over S and derive the expression for the eigenvalue probability distribution for non- 
Gaussian distributions of the form ( |1.1|) . In [Q this difficulty was circumvented by 
using the method of hermitization. 

As an explicit example, the case V{(j)''(j)) = 2m?(f)^ (() + g{(l)^ ipY was studied in detail 
in 0]. As should perhaps be expected in advance, the following behavior in the 
parameter space m^^g > was found: for positive, the eigenvalue distribution 
was disk-like (and non-uniform), generalizing Ginibre's work, but as = — was 
made more and more negative, a phase transition at the critical value 

occured, after which the disk fragmented into an annulus. The density of eigenvalues 
was calculated in |^ in detail. 

The paper is organized as follows: In Section 2 we summarize the "method of 
hermitization" We present (without derivation) the general algorithm for finding 



the density of eigenvalues associated with ( p..l[ ) which was developed in [Q, and also 
add some new insight into the mechanism behind the "Single Ring" theorem. We 
then formulate a novel simple criterion on the couplings in V{(j)''(f)) to decide whether 
the shape of the eigenvalue distribution is a disk or an annulus. Finally, we discuss 
some generic features of the disk-annulus phase transition. In particular, we prove 
that the Green's function associated with the hermitean matrix (which plays 
an important role in the "hermitization algorithm" just mentioned) is continuous 
through the disk-annulus phase transition. 

In Section 3 we provide simple new expressions for the outer radius Rout and for 
the inner radius (in the annular phase) of the eigenvalue distribution of the non- 
hermitean matrix 0, and for the corresponding boundary values p{Rout) and p{Rin) 
of its eigenvalue density, in terms of the moments 



< >= I daa'^pia) (k = 0, ±1, ■ ■ ■) 



j daa'^p{a) {k 

of the eigenvalue distribution p(cr) of the hermitean matrix 0^0. Thus, we find that 



and 



Rlut =<<y> 



P{Rout) ^ °" 



< > - < (T >2 ■ 



We see that -R^^^ is simply the average of cr, and the density p{Rout) is inversely 
proportional to the variance of a. 

Similarly, we find that in the annular phase. 



1 =(i 



and 

p{Rin) 



'in 



< a-^ > - < cr~^ >2 



Thus, Rj^^ is simply the a ^ moment of p{<y), and the density p{Rout) is inversely 
proportional to the variance of a~^. 
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In Section 4 we verify that the exphcit analytic expressions in Q concerning 
the quartic ensemble V^(0V) = 2m^0"''0 + g{(f)'^(f)Y are consistent with the results of 
Section 3. We also compare these analytic predictions with results of Monte-Carlo 
simulations of the quartic ensemble for various values of th? and g. The numerical 
results we obtained for the eigenvalue distribution in the disk phase and in the annular 
phase, as well as some quantitative features of the disk-annulus transition are in good 
agreement with the analytic predictions in [Q. 

The "Single Ring Theorem" may seem surprising at first sight. Our explana- 
tion of why the single ring theorem is not that surprising rests upon the simple 
argument that fragmentation of the eigenvalue distribution of into several dis- 
joint segments does not necessarily imply that the eigenvalues of trace out an- 
nuli obtained, loosely speaking, by revolving the segments of the eigenvalue distri- 
bution of 0^0 into the complex plane (see the discussion in Section 2). In Section 
5 we carry a numerical check of the "Single Ring" theorem for the sextic potential 
y(0V) = ?T2,^0"''0+|(0"'"0)^ + 3(0V)^) which is the potential of lowest degree for which 
the eigenvalues of cf^cj) may split into more than a single segment (in this case, two 
segments at the most). We generated numerically an ensemble in which the spectrum 
of (f^(j) is split into two separated segments, yet we found that the spectrum of is a 
disk, and not a configuration of a disk encircled by a concentric annulus, as one would 
perhaps naively expect by rotating the two-segment spectrum of 0V in the complex 
plane. 

In the Appendix we briefly review the multi-cut phase structure of matrix en- 
sembles with generic y(0^0), and then specialize to the phase structure of the sextic 
potential ensemble. 
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2 The Method of Hermitization and Non-Gaussian 
Ensembles 



Here we very briefly summarize the "method of hermitization" 0, Q in the form of an 
algorithm, followed by a general discussion of the phase structure of the eigenvalue 
distribution. 

Let us first introduce some notations and definitions. The averaged density of 
eigenvalues 

^) = E '^(^ - Re A,) 5{y - Im A,)) (2.1) 

i 

of the non-hermitean matrix 0, may be determined from the the Green's function 
associated with 0, namely 

G(^--)H>^>^/A'4^f^, (2.2) 

in terms of which[| 



p{x,y) = -d* G{z,z*). (2.3) 

TT 



The probability distributions ( |1 . 1|) studied in this paper are invariant under 
e*°0, rendering 

p{x,y) = p{r) /2tx (2.4) 

circularly invariant. Rotational invariance thus leads to a simpler form of the defining 
formula (p.2|) for G{z,z*) which reads 



7(r) = zG{z, z*) = J r'dr' p(r') , (2.5) 



whence 

p{r) = -^. (2.6) 
r ar 



^We use the following notational conventions: for z — x + iy we define d = — ^ — 
so that dz = 1. Similarly, we define = gfr = 5 + , so that d*z* = 1 and also d*{l/z) = 
■kS^^^z). Finally, we denote \z\ = r. 
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Clearly, the quantity 7(r), which can be thought of as the integrated eigenvalue 
density, is a positive monotonically increasing function, which satisfies the obvious 
"sum-rules" 

7(0) = and 7(00) = 1 . (2.7) 

In particular, observe that the first condition in ( |2.7| ) insures that no S{x)S{y) spike 
arises in p{x,y) when calculating it from (|2.3| ) with G{z,z*) given by (|2^ ), as it 
should be. 

It was shown in that by applying a simple trick, the desired Green's function 
of a non-hermitean random matrix could be obtained with a minimal amount of 
work, by judiciously exploiting the existing literature on random hermitean matrices. 
The algorithm, according to [^], for finding the Green's function and the averaged 
eigenvalue density of a non-hermitean random matrix drawn from a non- Gaussian 
ensemble P(0) = (l/Z)e-^t=-^(<^^'^) (Eq. (0)) is as follows: 

Start with the Green's function^ 

00 

F{w) = (-tr — -) ^ / i^^-L- , (2.8) 



'N '■^^ w — (p'(p' J w — a 
where 

P(/i) = ^(tV) 5(/x-0V)) (2.9) 

is the averaged eigenvalue density of 0^0. Then, the desired equation for 7(r) = 
zGiz^ z*) is 
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0. (2.10) 



Thus, given F one can solve for 7(r) using this master equation. 

Eq.( |2.10D is an algebraic equation for 7(r) and thus may have several r dependent 
solutions. In constructing the actual 7(r) one may have to match these solutions 
smoothly into a single function which increases monotonically from 7(0) = to 

^Of course, F(w) is already known in the literature on chiral and rectangular block random her- 
mitean matrices for the Gaussian distribution pi], p^ , [l^ , p^ , as well as for non-Gaussian probability 
distributions of the form (1.1) with an arbitrary polynomial potential F((/)'l'(/))[n8^, n9^, pO|. 
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'y{oo) = 1. An explicit non-trivial example of such a procedure is the construction of 
7(r) in the disk phase of the quartic ensemble^]. 

A remarkable property of ( |2.1(j| ) is that it has only two r-independent solutions: 
7 = and 7 = 10]. Since the actual 7(r) increases monotonically from 7(0) = to 
7(00) = 1, we immediately conclude from this observation that there can be no more 
than a single void in the eigenvalue distribution. Thus, in the class of models governed 
by -P(0) = ^e^'^^^^^'^'''^^ (Eq. ( [L.l|) ), the shape of the eigenvalue distribution is either 
a disk or an annulus, whatever polynomial the potential 1^(0^) is. This result is the 
"Single Ring Theorem" of g. 

The "Single Ring Theorem" may appear counter-intuitive at first sight. Indeed, 
consider a potential with several wells or minima. For deep enough wells, 

we expect the eigenvalues of to "fall into the wells". Thus, one might suppose 
that the eigenvalue distribution of to be bounded by a set of concentric circles of 
radii < ri < r2 < ■ ■ ■ < r^^^^, separating annular regions on which p(r) > from 
voids (annuli in which p(r) = 0.) A priori, it is natural to assume that the maximal 
number of such circular boundaries should grow with the degree of V, because V may 
then have many deep minima. Remarkably, however, according to the "Single Ring 
Theorem" the number of these boundaries is two at the most. 

To reconcile this conclusion with the a priori expectation just mentioned, note 
that while the eigenvalues of the hermitean matrix may split into several disjoint 
segments along the positive real axis, this does not necessarily constrain the eigenval- 
ues of itself to condense into annuli. Indeed, the hermitean matrix can always 
be diagonalized 0V = WA'^U by a unitary matrix U, with A^ = diag(A^, A^, ■ ■ ■ , A^), 
where the Aj are all real. This implies that = V'^AU, with V a unitary matrix as 
well. Thus, the complex eigenvalues of are given by the roots of det(2; — AW) = 0, 
with W = UVl Evidently, as W ranges over U{N) (which is what we expect to 
happen in the generic case), the eigenvalues of AW could be smeared (in the sense 
that they would not span narrow annuli around the circles \z\ = |Aj|.) 

The last argument in favor of the "Single Ring Theorem" clearly breaks down 
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when W fails to range over U{N), which occurs when the unitary matrices U and V 
are correlated. For example, may be such that W = UV'^ is block diagonal, with 
the upper diagonal block being a K x K unitary diagonal matrix diag(e*'^i, ■ ■ ■ , e*"^^"') 
(and with K a finite fraction of A^). In the extreme case = A^, in which W is 
completely diagonal, W = e'''^ = diag(e*'^i, ■ ■ ■ , e*"^^), we see that = We^'^AU is a 
normal matrix]^ with eigenvalues diag(e*'^^ Ai, ■ ■ ■ , c^'^^Aat). Thus, normal matrices, or 
partially normal matrices (i.e., the case K < N), evade the "Single Ring" theorem: if 
the first K eigenvalues , A^, ■ ■ ■ , A|f of 0V split into several disjoint segments along 
the positive real axis, the corresponding eigenvalues of will split into concentric 
annuli in the complex plane obtained by revolving those A-segments. Normal, or 
partially normal matrices are, of course, extremely rare in the ensembles of non- 
hermitean matrices studied in this paper, and do not affect the "Single Ring" behavior 
of the bulk of matrices in the ensemble. 

We end this section by showing how simple features of F{w) indicate whether 
the domain of the eigenvalue distribution is a disk or an annulus. As is well known 
ri| , pjj , for V a polynomial of degree p, the Green's function F{w) is given by0 



F{w) = -V'{w) - P{w)^{w - a){w - h) , (2.11) 

where 

p-2 

PM= ^kwK (2.12) 
fc=-i 

The real constants < a < 6 and Ck are then determined completely by the require- 
ment that F{w) ^ ^ as w tends to infinity, and by the condition that F{w) has at 
most an integrable singularity as — > 0. Thus, if a > 0, inevitably c_i = 0. However, 
if a = 0, then c_i will be determined by the asymptotic behavior for w large. 

According to the "Single Ring" theorem [Q, the eigenvalue distribution of is 
either a disk or an annulus. The behavior of F{w) as w ~ turns out to be an 
indicator as to which phase of the two the system is in, as we now show: 



^That is, [0, (f)^] = 0. 

"^Here we assume for simplicity that the eigenvalues of 0^ condense into a single segment [a, b] 
Discussion of condensation of 0^ eigenvalues into more segments appears in the Appendix. 
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A. Disk Phase: In the disk phase we expect that p(0) > 0, as in Ginibre's 
case. Thus, from (p.6|) p(r) = {1 / r) {dj / dr) = 2{d'j/dr^) and from the first sum rule 
7(0) = in ( ^.71 ) we conclude that 

7(r) ~ ^pioy (2.13) 
near r = 0. Therefore, for r small, (|2.1CI|) yields 

^(-^ + -)-^' (2.14) 

namely, F{w) ~ w ~ 0, as we could have anticipated from Ginibre's case. f\ 

This means that in the disk phase we must set a = in Consequently, in 

the disk phase c_i does not vanish. We can do even better: paying attention to the 
coefficients in and (|2.12|) (with a = 0) we immediately obtain from (|2.14|) that 



B. Annular Phase: In the annular phase 7(r) must clearly vanish identically in 
the inner void of the annulus. Thus, ( |2.10| ) implies that F{w) cannot have a pole 
at w = 0, and therefore from (p.ll| ) we must have c^i^/ab = 0. Thus, the annulus 
must arise for c_i = (the other possible solution a = 0, c_i 7^ leads to a disk 
configuration with 7 = only at r = 0, as we just discussed.) 

Thus, to summarize, in the disk phase F has the form 



Fdiskiw) = ^V'{w) - \J^w-^ + co + ciw + --- + Cp_2 wP-^j ^w{w-h) , 

(2.16) 

while in the annular phase it has the form 



F, 



annulus 



[w) = W\w) - (co + ci w + ■ ■ ■ + Cp_2 w^-^) ^{w-a){w-h) . (2.17) 



^In the Gaussian case, V = we have 2y/wF{'w) — ^fw — \Jw — 4, whence the roots of ( |2.10D 
are 7 = 0, 1 and . We note that 7 = is unphysical, 7 = (i.e., G — z*) corresponds to Ginibre's 
diskp2[, and 7 = 1 is the solution outside the disk. 
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Having determined F{w) in this way, i.e., having determined the various unknown 
parameters in (|2.16|) or in (|2.17|) , we substitute it into ( p.lO|) and find G{z,z*). We 
can thus calculate the density of eigenvalues p(r) explicitly for an arbitrary V. 

We now turn to the disk-annulus phase transition. An important feature of this 
transition is that F{w) is continuous through it. To see this we argue as follows: 
By tuning the couplings in V, we can induce a phase transition from the disk phase 
into the annular phase, or vice versa. Note, of course, that we can parametrize any 
point in the disk phase either by the set of couplings in V or by the set of parameters 
{c-i,co, ■ ■ ■ Cp-2',b} in ( p.l6|) . The "coordinate transformation" between these two 



sets of parameters is encoded in the asymptotic behavior of F{w). Similarly, we can 
parametrize any point in the annular phase either by the set of couplings in V or by 
the set of parameters {cq, Ci ■ ■ ■ Cp_2; a, b} in (|2.17| ). Due to the one-to-one relation (in 



a given phase, once we have established it is the stable one) between the couplings 
in V and the parameters in F{w) — ^V'{w) (namely, the c„'s and the locations of 
the branch points of F{w)), we can describe the disk-annulus transition in terms of 
the latter parameters (instead of the couplings in V). Clearly, the transition point is 
reached from the disk phase when p(0) = 0, that is, when c_i in (^.161) vanishes: 



„crit 



0. (2.18) 



Similarly, the transition point is reached from the annular phase when the lower 
branch point a in ( |2.17| ) vanishes. Thus, e.g., in a transition from the disk phase into 
the annular phase, F^iskiw) in ( p.l6|) would cross-over continuously into Fannuiusiw) 



in ( |2.17D through a critical form 

FUw) = lv:,uM - (C* + cr* W + --- + d^l^ w^-^) ^w{w - . (2.19) 

The continuity of F{w) through the transition was demonstrated explicitly in 0] 
for the quartic ensemble V{(f)^(f)) = 2'm'^(j)^(f) + g{(f)^(f)Y (see also Section 4). 

This discussion obviously generalizes to cases when F{w) has multiple cuts, which 
correspond to condensation of the eigenvalues of 0V i^ito many segments. If = 
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is a branch point of F[w), that is, if the lowest cut extends to the origin, we are in 
the disk phase. 



Fdisk{w) 



\v\w)-(c 



'^c-iw ^ + co + ci w -\ + Cp-2 ^) ^w{w - bi) ■ ■ ■ {w - hn) , 

(2.20) 

with < 61 < ■ ■ ■ < 6^1- The relation ( p.l5[ ) then generalizes to 



C-l 



P(0) 



(2.21) 



Since c_i must be real we conclude that such a configuration exists only for n odd. 
If the lowest branch point in F{w) is positive, we are in the annular phase with 

Fannulus{w) = (w) - (^Cq + Ci W ^ h Cp-2 W^^'^^ \J [w - a){w - hi) ■ ■ ■ {w ~ 6„) . 

(2.22) 

The phase transition would occur when the couplings in V{(f)^(j)) are tuned such that 
Fdiskiw) and Fannuiusiw) uiatch continuously, as was described in the previous para- 
graph. 
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3 Boundaries and Boundary Values 



Remarkably, with a minimal amount of effort, and based on the mere definition of 
F{w) (Eq. ( |2.8| ), which we repeat here for convenience) 

OO 

1 1 \ _ /" p{o')da 







we are able to derive simple expressions for the location of the boundaries of the eigen- 
value distribution and also for the boundary values of p{r) in terms of the moments 
of p(cr), which, we remind the reader, is the density of eigenvalues for a hermitean 
matrix problem. 



To this end it is useful to rewrite our master formula (|2.10|) for 7(r) as 

wF{w) = 7 (3.2) 

with 

w = -. 3.3 

7-1 

We start with the outer edge r = Rout (either in the disk phase or in the annular 
phase.) Near the outer edge 7 1 — , and thus w —00. We therefore expand F{w) 
in powers oil/w and obtain from (|3rT| ) - (|373|) 

i^, 1^2 

^3 



< (T > 7 - 1 n (7 - 1)^ , , 



where 



< a'^ >= / p{a)a''da (3.5) 



are the moments of p{a) (which is of course normalized to 1.) For the class of models 
we are interested in here, all the moments < cr'^ >, k > are clearly finite.^ Thus, at 
the outer edge r = Rout (where of course 'j{Rout) = 1), all terms with < a'' >, k > 2 
drop out of ( p.4|) and we obtain 

Rlut =<cr> . (3.6) 



^p(cr) = (l/7r)Im_F((T — ie) is supported along a finite segment (or segments), and its singularity 
at (T = is no worse than a^^^'^. 
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Namely, i^^Mt is simply the first moment of p{cr). 

We now calculate the boundary value p{Rout)- Approaching Rout from the inside, 
we substitute 7 = 1 — / and = R^outi.^ " ^) (with /, 5 << 1) in (P^. After some 
work we obtain / = ^^^^^^^2 ^ + 0{5'^), namely, 

7 = 1 Z""^' . 5 + 0{5^) . (3.7) 

Thus, from p(r) = 2{d^/dr^) (Eq. (|2;|)) and (g]!]) we find 



2^2 

P{Rout) = 5 — n ■ (3i 

< > _ < o- >2 



The density p{Rout) is inversely proportional to the variance of cr! 

For the p{a) under consideration here, < cr^ >, and consequently p{Rout), are 
always finite. Outside the boundary p(r) vanishes identically, of course, and thus, 
p(r) always "falls off a cliff" at the boundary, for all probability distributions of 
the form ( |1.1| ) with V polynomial. It would be thus interesting to study circularly 
invariant matrix ensembles P{(j)^(j)) such that the eigenvalue distribution p{a) of 
has a finite < a > but an infinite < o"^ >. Then p{Rout) would vanish. This would 
naturally raise the question whether in such situations, p(r) behaves universally near 
the edge (that is, if near the edge it vanishes like {Rout — with e being some 
universal exponent). We do not pursue this question further in this paper. 

We now turn to the annular phase, and focus on the inner edge r = R^ of the 
annulus. According to the discussion at the end of Section 2 (see Eq. (|2T7|) and the 



discussion above it), a > in ( p.ll| ), and thus F{w) is analytic in the domain |w| < a. 



Expanding (|3.1|) in powers of w, we obtain from 



-(-)=^(-^)+^'(-^) + ---- (3.9) 



— = w{ — )+W - 

A little above the inner radius, into the annulus, clearly 7 — and w 0— in ( p.3|) . 
Thus, setting w = in (|3.9| ) we obtain 

-Br = (-)- (3-10) 

Kin 
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R^^ is simply the moment of p(o"). 

We can now calculate the boundary value p{Rin)- Near the inner edge we parametrize 
= -Rf„(l + 5) with 5 « 1 (and of course, 7 << 1 to begin with.) Since p(cr) obvi- 
ously vanishes for cr < a, all moments < a^^ > in ( |3.9|) are finite. Thus, dropping all 
terms with < a"'^ >, A; > 3 in flS.QI), we obtain after some work 



7 = R^'—^-^—T^ + 0{6') . (3.11) 



< a-^ > - < o-^i >2 
It then follows from ( p. 61) and (|3.1CI|) that 

p(i?,„) = ^ . (3.12) 

The density p{Rout) is inversely proportional to the variance of cr~^ . 

From (|2.11|) (or ( p.l7| ) ) we learn that in the annular phase pannuiusi.^^) = ^Iiii-F(cr— 



it) = polynomial(cr) y (cr — a){b — a) , < a < a < b (and vanishes elsewhere.) Thus, 

b 

< >= S{p{a) / a^)da , k = 1,2 are finite. Therefore, Pannuius{o-) jumps from 

a 

zero (in the inner void of the annulus) to a finite value at the inner edge Rm- Note, 
however, that when a —>■ 0, that is, in the annular to disk transition, < 1/a > 
remains finite, but < l/o"^ > diverges like 1 / ^/a. (For a particular example see Eq. 
(|4.18|) .) Thus, from (|3.10| ) we see that Rinnner{o, = 0), the critical inner radius, is 



finite. The annulus starts up with a finite inner radius. Also, in this limit, we see 
from (|3.12| ) that p(-Rj„) vanishes hke ^/a. As we approach the annulus-disk transition, 
the discontinuity in p(r) at the (finite) inner edge disappears. 

We saw at the end of Section 2 (see Eqs. ( p.l6| )-( pTT9D ) that F{w) is continuous 
through the disk-annulus phase transition. Thus, our master formula wF{w) = 7 to 
determine 7(r) (Eq. (|3.2|) ) is also continuous through the transition. Consequently, 
p(r) = {l/r){d'-f /dr) must remain continuous through the disk-annulus transition, and 
has (at the transition) the universal behavior described in the previous paragraph. 
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4 Phase Transitions in the Quartic Ensemble 



The disk-annulus transition in the quartic ensemble 

= 2mVV + 9{<P'^<PY (4.1) 

was studied in detail in 0]. The annular eigenvalue distribution Pannuiarij) and the 
disk eigenvalue distribution pdisk{r) for this ensemble were calculated explicitly in 
According to the expressions given in 0], as the critical point is approached from the 
annular phase, Pannuiar{r) behaves precisely as described in the paragraph following 
Eq. (|3.12|) at the end of the previous section (see also Eq. ( |4.13| ) below, at /i = pc-) 
Also according to 0], as the critical point is approached from the disk phase, Pdiskij') 
gets completely depleted inside a region of radius Rin{pc) (remaining continuous at 



r = RiniPc)- See Eq. ( [4.20|) below.) Thus, p(r) for the quartic ensemble is continuous 
through the disk-annulus transition. 

In this section we verify the expressions (|3^), ( |3.8| ), ( p.lO| ) and (|3.12|) for Rout, 



p{Rout), Rin and p{Rin) for the quartic ensemble ( [4.1|) against the explicit expressions 
for these quantities given in and also provide ample numerical results concerning 
the disk phase, the annular phase, and the transition between them, in support of 
the analytical results. In what follows we have omitted many technical details that 
can be found in [Q. 

4.1 The Disk Phase 

For m? > —\f2g the density of eigenvalues is a disk. According to we have 

F{w) = + gw - + ^Jw{w - b) (4.2) 

with 

2m2 + y/m* + 6g , , -2m^ + 2y/m* + 6g 

c = , and b = . (4.3) 

3 3g 

According to Eqs. (5.8) and (5.9) in [Q, the eigenvalue density in this phase is 

/ N ^ 2 . 2 bc^ - (m^ + 2gr^)\l + 2(m^r^ + gr^)] ,,,, 
Pd^skir) = 2m^ + 4gP + 2[sgn (- - r^)] / ^ ^ ^ ^ (4.4) 

4 ^'[l + 2(m2r2 + gr^)Y - Abc^r^ 
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inside a disk of radius Rout, where 

2 bc^ — 2m^ (m*^ + 6(7)^/^ — m^im^ + 9(7) 



Thus, from ( [4.4| ) and ( |4.5| ) we have [] 

P'-'^''"-'^ 2g - bc\bc^ - 2m2) 1 - 2RUgRl^, + m^) " ^^"'^ 



These resuhs should be compared with the predictions of Section 3. From ( ^.2| ) 
we can read-off the density of eigenvalues p{a) = (l/7r)ImF((T — ie) of (p'^cp as 



1 /c 



p(a) = - - + ^ ^aib - a) (4.7) 

for < 0" < 6, and zero elsewhere. We can readily check that ( |4.7| ) is properly 
normalized to 1. 

The first two moments of ( ^4.7| ) are 

^_ 1 9b\ _ {m' + Qgf/'-m\m^ + 9g) 

2 1 2 i r 2 i 27(72 ' ^ ■ 



and 

4 



2 lfb\ fSc \ 27g^ + 18gm^ + 2m^ -2m'^{6g + my/^ 



Thus, 



< o-^ > - < 0- 



63 



46c^ + bg{b'^g - 10) + 4c(6^^ - 4) 



256 

297^3 + 108^W - IS^m^ - 4mi2 - 2m2(9^ - 2m'^)(6^ + m^fl'^ 



. (4.10) 



1458C/4 

Comparing ([4.5|) and (^78|) we immediately verify (|3.6| ), -R^ut =< o" >• After some 
additional work, using ( |4.10| ) and ( [4.5| ) in (|3l8|), we can see that (p.8|), namely, that 
p{Rout) = 2-Ro„i/(< > — < a >2), coincides with 



''It is straight forward to show that sgn (| — ^ottt) — ~1 ([4. 4]): just substitute m** = A^cos^ 



^4 

and 6g = A^sin^a and obtain that Rl^^ -\ = {A^ /l()ig^){l - cosaf 
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4.1.1 Numerical Results for the Disk Dhase 

We have generated numerically random matrix ensembles corresponding to the quartic 
potential (|4.1| ) in the disk phase, for = 1 fixed and for various values of the coupling 
g (and for various sizes of matrices), and measured pdiski'f') for these realizations. 

The generation of the matrices was done by a standard Metropolis Monte Carlo 
approach. A random change had been suggested in the real and imaginary parts of 
one of the elements of and then the change in V{(j)) was evaluated. This "move" 
was accepted unconditionally if V were decreased, and with probability p = e"^^ if V 
were increased. General theorems on Monte Carlo then guaranteed that the resulting 
probability distribution of was the desired one. After the matrices were generated, 
their eigenvalues were determined with a standard solver from the LAPACK library. 
We tuned the size of the suggested changes in (p so that the acceptance rate was 
about one-half, and monitored the equilibration and autocorrelation times to ensure 
our starting configuration had evolved properly and error bars have were accurate. 
In particular, the local changes in (p made the matrices correlated over some number 
of random changes, however, local changes also allowed one to employ various tricks 
to evaluate the change in V rapidly. 

In Figure 1 we display our numerical results for pdiski'f') for 128x128 dimensional 
matrices, and compare them to the analytical large- iV result ( [4.4|) of . (As a trivial 
check of our numerical code, we also included in this figure the results for the gaussian 
(Ginibre) ensemble.) 

Evidently, the agreement between the numerical and the analytical results is good. 
Note the finite- effects near the edge of the disk. 

4.2 The Annular Phase 

For < — the stable eigenvalue distribution is annular. For convenience, let 
us switch notations according to = — /i^, and also write /i^ = y/2g. 
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According to 0] we have 

F{w) = + gw — g \J {w — a){w — h) (4-11) 

with 

a=^-,/?, and b = ^ + . (4.12) 
9^9 9^9 

We see that a = (2//x^)(/i^ — /i^) which is positive for /x^ > /x^, as it should be, by 
definition. 

According to Eqs. (5.16) - (5.19) in [^], the eigenvalue density in this phase is 

Pannuiusir) = 8g (^T^ '^^=^9 (t' - (4.13) 
inside an annulus Rm < r < Rout, where 



/i2 + v//x4 _ 2g _ f^' + y/i^ - /i^ 
2y " l^t 



RL = " " = (4.14) 

and 



Rlt = T = ^ • (4.15) 

9 H'c 



Thus, we see immediately that 



and 

piRout) = 4/i^ (4.16) 

Note that p{Rin) = at /i = fic, as expected. Also note that the critical annulus has 
a finite inner radius: Rf^^pc) = ^/pl > 0- 

We now compare these results with the predictions of Section 3. From ( [4.11 
read-off the density of eigenvalues of 



we 



~pia) = ^-^{a-a){h-a) (4.17) 



for a < a < 6, and zero elsewhere. We can check that (f4.17|) is properly normalized 
to 1. 
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The relevant moments of (|4.17|) are 



< cr > 



< a" > 



g f a + h\ fh — a\ fi"^ 
2 I ^ J I ^ J = 7 



g fh — a" 
2 



' h + a \ 1 
h-a] ^4 



< - > 
a 

and 



< — > 



-iVh-V^f = ^?-J^^-2g 



2 ^ - 2(7 



(4.18) 



Comparing ( |4.14| ), ( ^.151 ) and the first and third equations in ( [4.1^ ), we verify 
(|3.6|) and ( p.lO|) straightforwardly. 



Further, we find from ( 4.18 ) that 



2 2 
< O > - < O > = — 



and 



1 



< ^ > - < - > 
Thus, comparing with ( [4.16|) we find that 



2^7 

g-2fi'' + fi 



(4.19) 



< > _ < o- >2 



4//^ = p{Rout) 



and 



2i?„ 



< > - < >2 



4 J/i^ - /i^ = p(i?,„) 



and verify ( |3.8| ) and ( |3.12D for the annular phase. 



4.2.1 Numerical Results for the Annular Phase 

In Figures 2.a-2.c we display our numerical results for Pannuiusir) for matrices of 
various sizes, and compare them to the analytical large- iV result ( ^.13| ) of Q]. In 

these figures we hold /i^ = 0.5 fixed, and increase g from 0.025 to 0.10. 

®Here we have /i^ = 0.5 = ^^^/2^/2g. Thus increasmg g as indicated in the text brings us closer 
to the disk-annulus phase transition. 
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4.3 The Disk-Annulus Phase Transition 

The phase boundary separating the disk phase and the annular phase in the — g 
plane is the curve = —\/2g. 

Consider approaching this boundary from within the disk phase by setting rn? = 
— ^/2g + 6, with 6 positive and small. Then, using (|4.3| ), we find to first order in 6 that 



c = S/2 and b = 2y(2/g) — 5/g. In particular, at the phase boundary itself c = 0, 
in accordance with ( p.l8|) and ( p.l9|) . It was shown in 0] that as one approaches 
the critical point = —\f2g from the disk phase, the density of eigenvalues of 
approches the particularly simple critical configuration 

r , r2<l/^ 

Pcrit(r) = <^ (4.20) 

Thus, as we decrease 5 to zero, pdisk{r) (Eq. ( W^ ) becomes increasingly depleted 
inside the disk < b{6)/4, reaching complete depletion at 5 = 0, at which point the 
disk breaks into an annulus. We also note that at the phase boundary ([4.2|) reads 



F{w) = -J2g + gw - g 



\ 



w(w- 2^\ . (4.21) 



Consider now approaching the phase boundary = —\/2g from within the 
annular phase. Thus, we set /i^ = \f2g + 5, with b positive and small. Then, since all 
the expressions in ( [4.12|) are linear in /x^, we find that a = 5 / g and h = 2^ {2/ g) + 5 /g. 
In particular, at the phase boundary itself a = 0, and h = l/y/2g. Therefore, at the 
phase boundary ( |4.11| ) reads 



Fi^) = -^'^9 + - g 



( /2' 
A w \ w — 2\ — 



which coincides with ( [4.21|) . Thus, F{w) (and consequently, the eigenvalue density 



of 0^) is also continuous at the transition, in accordance with ( |2.19| ). 

Note from ([4.14|) and ([4.15|) that at the transition R'^^^ = l/fi^ = 1/ ^/2g is finite, 



and coincides with the radius (squared) of the depleted region in the disk configuration 
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( |4.20| ). Right at the transition, the disk breaks into an annulus with a finite hole! 



Note also that R^^^ = 2//i^ = y^/g, which coincides with the disk's R^^^ at the phase 
boundary. Thus, at the phase boundary /i^ = /i^ ( 4.1^ ) coincides with ( [4.2U| ), namely. 



p(r) is continuous at the transition from the disk phase to the annular phase. 
4.3.1 Numerical Simulation of the Disk-Annulus Phase Transition 

We have measured the density of eigenvalues p{r) of matrices of size 128x128, taken 
from the the quartic ensemble with /i^ = 0.5 and for g = 0.025,0.05,0.1,0.125,0.15 
and 0.175. The results are displayed on Figure 3. 

For these values of g, we start in the annular phase at the lowest value of g. For 
our set of parameters we have /i^ = 0.5 = [x^^j^^/^g. Thus, increasing g (while keeping 
/i^ fixed at 0.5) brings us closer to the disk-annulus phase transition, which occurs 
(at large A^) at g^ = 0.125. Increasing g beyond that, puts us into the disk phase. 

The first three profiles on the right in Figure 3 belong to the annular phase. Their 
behavior is consistent with our discussion in Section 4.2 of the annular phase. Indeed, 
as g increases towards the transition point at gc = 0.125, these three graphs exhibit 



the expected decrease of Rf^ = (/i^ + y/i'^ — lit)/ iJi^ (Eq. ( [4. 141) , with /x^ = y/2g) and 
the decrease of Rl^^ = 2/iV/i^ (Eq. (|4l5|) ). 

The critical density profile, corresponding to gc = 0.125, is the fourth profile (from 
the right). For our choice of parameters, the theoretical boundary radii of the critical 
annulus, i.e., at = 0.125, are i?^"* = 1/fic = and i?^"/ = V^/f^c = 2. These 
boundary values fit nicely with the features of the critical profile in Figure 3. 

Finally, the last two profiles in Figure 3 have pronounced tails extending to r = 
and thus belong to the disk phase. 
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5 Phase Transitions in the Sixth Order Potential 
and the "Single Ring" Theorem 

The sextic potential 

is the potential of lowest degree in ( |1 . 1| ) for which the eigenvalues of 0V niay split 
into more then a single segment. In fact, it is easy to see that there can be at most 
two eigenvalue segments in the spectrum of 0^0. 

The qualitative features of the support of the eigenvalue density associated with 
(|5.1|) can be deduced by moving the cubic V{x) = rri^x + + around in the 
plane (i.e., by varying its couplings, fixing, say g = 1), and concentrating on x > 0. It 
is obvious from such considerations that there are three qualitatively different phases 
in the spectrum of 0^0. In two of the phases the eigenvalues Sj of live in a 
single segment. In one of these single segment phases, the segment includes the origin 
(0 < s < a), but in the other it does not (0 < a < s < 6). Following the general 
discussion in the last part of Section 2, we would expect that the spectrum of (p itself 
is a disk in the first case, and an annulus in the second case. 

In the third phase of 0^0, there are two segments, one of which hits the origin 
({0 < s < a} \J{b < s < c}). (There is no two-segment phase of 0V which does not 
include the origin.) Thus, according to the discussion in the last part of Section 2, the 
non-hermitean matrix is expected to be in the disk phase in this case (rather than 
having its eigenvalue fill in a disk surrounded by a concentric annulus), in accordance 
with the "Single Ring" theorem. 

In this short section we limit our discussion to the two-segment phase of (A 
rather detailed sketch of the analytical conditions that determine the whole phase 
structure of the sextic ensemble ( ^.1|) is given in the Appendix.) Our purpose here is 
to demonstrate numerically the "Single Ring" Theorem for the eigenvalue distribution 
of matrices taken from the sextic ensemble ( |5.1|) . To this end, we have to identify 
points well within the phase in which the eigenvalues of 0V split into two segments. 
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We used the formalism of the Appendix to choose two ensembles in the two- 
segment phase of 0^0, for which we verified that the eigenvalues of formed a disk. 
The results for these ensembles are displayed in Figures 4 and 5. 

Figure 4 shows the scatter plot of the eigenvalues of together with the density of 
eigenvalues p(s) of 0^0 for ( p.l| ) with couplings = 7.372, A = —6.116 and g = 1.372. 
As can be seen on the right part of Figure 4, for these couplings, the eigenvalues of 
0^0 live in two separated segments: {0 < s < a = 1} \J{b = 2 < s < c = 3}. The 
solid line there is the large N theoretical curve, which was plotted according to the 
analysis we have described in the Appendix (see the discussion following (|A.2(]|) ). 
Evidently, the spectrum of is a disk, despite the spht support of p(s), in accordance 
with the "Single Ring" Theorem. 

Figure 5 is similar to Figure 4, but for (|5.1|) with couplings = 6.403, A = 4.184 
and g = 0.713, for which the eigenvalues of 0^0 live in the segments {0 < s < a = 
1} U{^ = 3 < s < c = 4}. The spectrum of remains a disk, even though the two 
segments of the support of p{s) are more separated than in Figure 4. 
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Appendix A The Phase Diagram of a Generic V{(f)^(f)) 



In this Appendix we briefly review the necessary theoretical aspects of multi-cut 
phases of 0^0. The first part of our discussion will apply for a generic polynomial 
V{(f)''(j)). Then, in the second part of the Appendix, we will specialize to the sextic 
potential 

For practical reasons, we eliminate some (or all) of the couplings in V in 
terms of the end-points of the segments containing the eigenvalue distribution of 
and use the latter as (part of) the coordinates in the phase diagram. In this way we 
can find rather easily which couplings in V{(f)''(j)) are needed to generate an eigenvalue 
distribution of (p^cp with a prescribed set of support segments. 



A.l A generic potential 



ReF{s - ie) = -V'{s) . 

2 



By definition (see Eq. (Ol 



1 1 1^1 

(where Si are the eigenvalues of 0V)- Thus, as usual, 

1 JL. 1 



F{s - te) = P.P^^( ) + ^vrp(s) 

iV j^-^ S Sj 



The saddle-point equation governing the Dyson gas of eigenvalues of (p'cj) is [19, 20 



(A.l) 



(A.2) 



(A.3) 



where p(s) is the density of eigenvalues of 

In order to study multi-cut configurations of p{s), we also need the auxiliary 



function [21 



G{s)= J dfi{V'{fi) -2F{fi-ie)) . 

ai>0 



(A.4) 
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In the last equation ai is the lowest branch point of F{w). Thus, from ( |A.3| ) and 
(|A.1|) , we see that for s real and in the support of eigenvalues, 

s 

G{s) = -2m J p(/i)rf/i (A.5) 

11 

is pure imaginary. — Im G{s) = 2Tip{s) is then positive and monotonically increasing 
(and reaches 27r when s hits the largest branch point). 

Stability of multi-cut distributions How do we know that a given distribution 
of eigenvalues is stable against migration of eigenvalues from one place to another? 
To answer this question, consider the Dyson gas energy functional 

Seff[p]= J p{s)V{s)ds-^ J p{s)p{p)\og {s - pY dsdp . (A. 6) 

s>0 s,/i>0 

A general variation of ( |A.6| ) under p(s) p{s) + 5p{s) is 

5Seff[p\= j V{s)6p{s)ds— J p{p)\og {s — p)'^6p{s) dsdp . (A. 7) 

Moving an eigenvalue from Sj to sj corresponds to Sp{s) = {l/N)[S{s — Sf) —S{s — Si)]. 
Thus, from ( [A.7| ) and (|A.4| ) (and after some work) we can show that such a move 
costs 

ASeff[p] = j^[G{sj)-G{s,)] (A.8) 



in energy [^. Such a rearrangement of eigenvalues costs energy only if 

Re A5e//[p] >0, (A.9) 



and therefore (|A.9| ) is the stability condition against such a rearrangement. Thus, a 
multi-cut F{w), where the eigenvalues coalesce into n segments 

[ai, 02] |J[a3, 04] U ■ ■ ■ lj[«2n-l, a2n], 

would be stable against migration of eigenvalues between neighboring cuts if and only 
if ( [A.9[ ) would hold for all neighboring pairs of cuts, and in both directions. Since 
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G{s) is real on the segments on the real axis that connect the cuts, this stability 
condition means 



G(a3) = G(a2), G{a,) = G{a^), ■ ■ ■ G{a2n-i) = G{a2n-2) ■ (A.IO) 
In addition, of course. Re G{s) < cannot happen anywhere for s > 0. The n — 



1 equations (|A.10|) comprise the desired stability condition for such an eigenvalue 
distribution. In addition to these conditions, we have to make sure that along the 
cuts themselves — Im G{s) > 0, which is just the condition that p{s) be positive. 

The n — 1 equations ( |A.1U| ), together with the obvious analytic properties of F{w) 
and its asymptotic behavior 

F(w) ~ - (A.ll) 
w 



as w ^ oo, determine F{w) uniquely. Indeed, as is well known, for a generic 
in view of ( |A.1| ) and ( [A.ll| ) (and as we discussed at the end of Section 2), F{w) (with 
n cuts) must be of the form 



F{w) = -V'{w) - P{w) 



2n 



\[(w-a,), (A.12) 



2 ' ' ' '%\ 

where 

dcg V—n—l 

Here ai < 02 ■ ■ ■ < a2n and c_i 7^ only if ai = (see Section 2). If ai > 0, then 
c_i = 0, and thus in such a case, there are 2n{a's) + {degV — n){c's) = n + degV^ 
independent parameters in the expression (|A.12|) for F{w). On the other hand, there 



are deg + 1 conditions from the asymptotic behavior ( [A.ll|) plus additional n — 1 
conditions from ( [A.10|) , which comprise a total of deg\^ + n conditions, equal to the 



number of unknown paramters. This balance remains if c_i appears in the game as 
an unknown parameter, because then oi = 0, so that the number of parameters does 
not change. Finally, we have to remember to impose the positivity constraint 

p{x) = -Im F{x -ie) >0, (A.14) 

71 
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which translates into a set of inequahties among the a's and c's. 

A convenient local parametrization of the phase diagram Recall, that the phases 
of (j)^(f) are specified by the number of segments in the support of p(s), i.e., the number 
of cuts in F{w) (and whether these cuts have w = as a branch point or not.) 
Thus, instead of the usual description of the phase structure in terms of the degV^ 
couplings in ^^(0^0), our strategy is to use deg\^ parameters out of the 2n branch- 
points ai, ■ ■ ■ , a2n of F{w) and the degP coefficients Ck (with the total number degV^ 
first saturated by the a's in ascending order), which we refer to as ^^phase coordinates" , 
to express (in a given phase) the couplings appearing in y(0"''0) (such as m?,\ and 
g in (|5.1|)), as well the as the c^'s and a^'s complementary to the phase coordinate 
parameter set. (See our discussion of the sextic potential below for concrete examples 
of this paramtrization.) 

We have to be careful in giving the expressions for, say, the couplings of V , in 
terms of the phase coordinates. This because for a given configuration of F{w), the 
equations from which we are to eliminate the couplings of V (such as the triad m^, A 
and g in ( ^.1|) ) as functions of the phase coordinates may have several solutions (in 
other words, the couplings in V are generally multivalued functions of the phase 
coordinates in a given phase). Thus, in a given phase, we must of course choose the 
parametrization of couplings in V which yields the minimal Seff[p] appropriate for 
that phase. 

This alternative parametrization is more convenient for our purposes in Section 
5. Indeed, once we are successful in expressing the couplings in V as functions of 
the phase coordinates, it will be very easy for us to tune the couplings in V{(j)^(j)) 
to a generic point in a given phase and also to approach the phase boundaries in a 
controlled manner. In particular, phase transitions appear here, for example, when 
branch points collide and become equal (at some common real positive value a). This 
process removes 2 a's and thus closes one cut {n —>■ n — 1) but adds an additional 
term to P{w). The number of unknown parameters drops by 1, but so does the 
number of stability conditions ([A.IOI). In the other dircetion, we can obviously reach 
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the same coexistence point, by tuning the parameters of P{w) to a point where it 



develops a hnear factor {w — a) = y (w — a)^ (with a > 0). Obviously, when these 
alternative phase coordinates approach a point on the coexistence surface from two 
different sides of the pahse transition, the respective sets of couplings of V, expressed 
as sets of functions of the two phase coordinate patches, coincide. Thus, they lead 
to the same Seff[p], which means that such a point is indeed a point on the phase 
boundary. 

A. 2 Results for the sextic potential 

From (|A.12 ) and ( |5.1| ), the general form of F{w) is 



F{w) = - ym^ + Xw + gw'^j — P(w)y polynomial . (A.15) 

A. single cut, disk phase: Here 

F{w) = i (m"^ + Xw + gw^) - + ^ + \Jw{w - a) . (A.16) 

There is a single cut, so ( |A.10| ) is trivial in this case, and ( |A.14| ) holds manifestly. We 
need only impose ( |A.11| ). In the end, we find 

16 - 8as - 2aH 



u 



a3 



g = 2u 

, ^ AaH + 8as - 16 

A = 2t — au = and 

9 a^u 8as — aH — 8 / , , „n 

m^ = 2s-at = . (A.17) 

4 2a ^ ' 



The phase coordinates are a, t and s. 

B. single cut, annular phase: We have 

F{w) = - (m^ + Xw + gw"^^ - {t + uw) sj{w - a){w - h) . (A.18) 



Here < a < 6. Again, there is a single cut, so ( [A.1(J| ) is trivial in this case too, and 



also ( [A.14| ) holds manifestly. We need only impose ( [A.ll|) . In the end, we find 
16 - 2t(a - 6)2 



u 



(a + 6)(a-6)2 
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g = 2u 

X = 2t-u{a + b)= At 

{a - by 
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(a - by 



and 



m 



u — t{a + b) 



a + b 



. + Qab + 6^ 
2(a + b) 



t- 



(A.19) 



The phase coordinates are a, b and t. 

C. two cuts, disk phase: We have 



F{w) = - (m^ + \w + gw'^^ — ( h t \ \Jw{w — a){w — b){w — c) . 



(A.20) 



Here < a < 6 < c. Note that in this case we can trade the three couphngs m?, A 
and g for the three branch points a, b and c. In this case there are two cuts, so for 
the first time (|A.10|) is not trivial. We first impose (|A.11| ). We find 

8 



~ a2 + 62 + c2 - 2{ab + ac + be) 
la!^ + b'^ + c" - a^{b + c) - b^{a + c) - c^{a + b) + 2abc 
2 a2 + 62 + c2 - 2{ab + ac + be) 

g = 2t 

X = 2s — t{a + b + c) and 

n / , N + 6^ + — 2(a6 + ac + 6c) 
= -(a + 6 + c)s -t . 



t 



(A.21) 



We have yet to impose ( |A.10 ), which is why t was not ehminated yet. Before doing 
that, we impose ( |A.14| ). Our conventions are always to take each cut from the appro- 
priate branch point to the left on the real axis. Thus, after some work, we find from 



TTpi^x) 



Im F{x — it) 



(- + t)Jx{a — x){b — x){c — x) < X < a 



+ (f + t)\/xix -a){x- b){c - x) b<x<c (A.22) 
otherwise . 
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We have to impose ([A.14| ) on ( [A.22|) . This means 



+ t<0 for 0<x<a 



Thus, we must have^ 



X 



- + t > 

X 



t > and 



and 

for b < X < c. 



bt < s < -at < , 



(A.23) 



where s is given in ( |A.21| ) . We are now ready to impose ([A.10| ) . Here it simply means 
G{a) = G{b), namely, 

b 

'(- + t)^Jx{x - a){b - x){c - x) dx = . {AM) 



Note from ( |A.23| ) that —b<s/t< —a, and thus the factor multiplying the square 
root in ( [A.24| ) flips its sign in the integration domain, so that the integral on the LHS 
of ( |A.24| ) may vanish. The latter equation may be expressed in terms of the elliptic 
integrals 

b 



I{a, b, c) 



^Jx{x — a){x — b){x — c)dx 



and 

J(a, b,c)= I \Jx{x — a){x — b){x — c) — . 



X 



(A.25) 



Following the usual procedure, we may express / and J in terms of complete elliptic 
integrals in a straightforward manner. (We do not bother to write these expressions 
here, since for our purposes in Section 5 we evaluated I and J numerically.) 
Finally, substituting these expressions in ( |A.24] ) we obtain 



s{t, a, b, c)J{a, b, c) + tl{a, b,c) = 



(A.26) 



^It is straightforward to check that the foUowing inequahties hold for the ensembles corresponding 
to Figures 4 and 5 in Section 5. 
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which we solve for t (recall from ( |A.21| ) that s is merely linear in t, and also that 
t > 0, in view of (|A.23|) .) Once t{a,b,c) is known, we can go back to ( |A.21| ) and 
evaluate m^(a, 6, c), A(a, 6, c) and g{a,b,c) explicitly. Our phase coordinates in this 
case are thus a, b and c. 
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N=128 
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Figure 1: Comparison between Monte-Carlo measurements of the density of eigen- 
values p(r) of matrices (f) of size 128x128, taken from the quartic ensemble V{(f)^(f)) = 
2m?(j)^(f) + g{(j)^(j))'^ with = 0.5 (disk phase) and for g = 0,0.5, 1,2,4 {g increases 
from bottom to top), compared to the analytical results of 0] (solid lines). At g = 
we obtain Ginibre's Gaussian ensemble with V = with its unit disk of eigenval- 
ues. 
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N=64,128,256 
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Figure 2: Results of Monte-Carlo measurements of the density of eigenvalues p(r) of 
matrices of sizes corresponding to = 64, 128 and 256, taken from the quartic 
ensemble with = —fi"^ = —0.5 (annular phase) for various values of the quartic 
coupling: g = 0.025 in [a], g = 0.05 in [b] and g = 0.1 in [c]. These are compared 
to the analytical results of 0] (solid lines). As N increases, the numerical results 
converge monotonically to the analytical results. 
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N=128 




Figure 3: Monte-Carlo measurements of the density of eigenvalues p{r) of matrices 
(f) of size 128x128, taken from the the quartic ensemble with /x^ = 0.5 and for g ~ 
0.025,0.05,0.1,0.125,0.15 and 0.175 {g increases from right to left). The first three 
profiles on the right (corresponding to the three lowest values of g) evidently belong 
to the annular phase. The fourth density profile from the right is the critical one 
(corresponding to Qc — 0.125). Finally, The last two profiles (which correspond to 
the two higher values of g) belong to the disk configuration. 
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Figure 4: Scatter plot of the eigenvalues of matrices of size 32x32, taken from 
(|5.1|) with = 7.372, A = —6.116 and g = 1.372 (left), and the corresponding 
density of eigenvalues p{s) of 0V (right). The solid line on the right is the analytical 
curve corresponding to ( |A.22 ). The support of p{s) is split into the two segments 
{0 < s < a = 1}U{^ = 2 < s < c = 3}, while the support of p(r) on the left is 
manifestly a disk. 
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Figure 5: Similarly to Figure 4, but with = 6.403, A = 4.184 and g = 0.713. The 
support of p{s) is split into the two segments {0 < s < a = 1} \J{b = 3<s<c = 4}, 
while the support of p(r) on the left remains a disk. 
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